Extrinsic Curvature, Geometric Optics, and Lamellar Order on Curved Substrates 
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When thermal energies are weak, two dimensional lamellar structures confined on a curved sub- 
strate display complex patterns arising from the competition between layer bending and compression 
in the presence of geometric constraints. We present broad design principles to engineer the geome- 
try of the underlying substrate so that a desired lamellar pattern can be obtained by self-assembly. 
Two distinct physical effects are identified as key factors that contribute to the interaction between 
the shape of the underlying surface and the resulting lamellar morphology. The first is a local 
ordering field for the direction of each individual layer which tends to minimize its curvature with 
respect to the three-dimensional embedding. The second is a non-local effect controlled by the 
intrinsic geometry of the surface that forces the normals to the (nearly incompressible) layers to 
lie on geodesies, leading to caustic formation as in optics. As a result, different surface morpholo- 
gies with predominantly positive or negative Gaussian curvature can act as converging or diverging 
lenses respectively. By combining these ingredients, as one would with different optical elements, 
complex lamellar morphologies can be obtained. This smectic optometry enables the manipulation 
of lamellar configurations for the design of novel materials. 

PACS numbers: 61.30.-v, 61.30.Hn, 02.40.-k 
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I. INTRODUCTION 



Though many are taught that there are merely three 
phases of matter, solid, liquid, and gas, the under- 
standing of broken symmetries, Nambu-Goldstone modes 
[U, 0, and the Anderson-Higgs-Kibble 's', l", '5^ mecha- 
nism allows prediction, control, and elucidation of other 
novel forms of matter. Liquid crystalline phases interpo- 
late between the simple forms of matter sketched above. 
The nematic has partially broken rotational invariance 
, the smectic phase has broken one-dimensional trans- 
lation invariance Q, and the hexatic phase Q has bro- 
ken two-dimensional rotational invariance. These mate- 
rials not only afford a deeper understanding of condensed 
matter, they (and their lyotropic cousins, diblock copoly- 
mers Q), are also of great technological interest, from 
displays to coatings, from drug delivery to hybrid ma- 
terials. Technological applications provide challenges for 
the theorist and, in particular, pose problems with imper- 
fect boundary conditions and geometries. In this paper, 
we focus on striped or smectic phases on curved sub- 
strates. Our work is motivated by elegant experiments 
by Hexemer and Kramer [To| who probed films of diblock 
copolymers on corrugated surfaces. Explaining the prop- 
erties of columnar layers in this geometry is important 
for controlling the resulting microstructure of the phase 
on frozen, undulating surfaces. 

Lamellar, smectic, striped, and columnar phases can 
all be treated in a single unified, framework on a two- 
dimensional surface, with the columns, layers, etc. lying 
in the tangent plane. The physics of smectics is beautiful 
and intricate: the coupling between their geometry and 



rotation invariance ensures the presence of essential non- 
linearities in the strain that lead to anomalous elasticity 
pll . [T^ and dramatic departures from linear elasticity 
even for small strains [l^, [3l . These nonlinearities make 
the theory of the nematic-to-smectic-A transition notori- 
ously difficult as well 15] and the effort to capture them 
has enhanced our understanding of global versus local 
symmetries. 

Posed with complex boundary conditions, smectics of- 
ten form focal conic domains [Toj in which the layers 
remain equally spaced but acquire a large curvature to 
accommodate their growth in a confined geometry. In 
the past decade or so, smectic phases with cubic order 
have been discovered [l7| , which have pushed our under- 
standing of the competing tendencies of uniform layer 
spacing and non-vanishing curvature. Indeed, the two 
requirements cannot generically be reconciled with the 
layer topology [l^. This leads to a subtle type of geo- 
metric frustration in which there is no local obstruction 
to finding uniformly-spaced layers but their construction 
leads to diverging layer curvatures. These focal conies 
are, in fact, analogous to the caustic singularities in op- 
tics [l^. It has been conjectured that the frustration 
in smectics can be lifted by introducing curvature in the 
background space |18| , much as the nematic blue phase is 
not frustrated on the surface of a four-dimensional sphere 
[2C|. It is, of course, not possible to experimentally study 
layered systems in high dimensional space. However, a 
simpler but not less subtle form of geometric frustration 
also exists in two-dimensional smectics lying on curved 
substrates [13, HH, [12] ■ In contrast to the case of two- 
dimensional crystals, [1^ HH the connection between 
strain and substrate curvature is indirect, and conse- 
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quently, configurations with uniformly-spaced layers are 
numerous. However, the layers nevertheless inherit cur- 
vature from the underlying substrate due to compression 
elasticity. Smectics exhibit quite different phenomena 
from other ordered phases on curved surfaces and new 
theoretical ideas are required to understand them. 

In this paper, we study configurations of uniformly- 
spaced layers on a curved surface, elucidating and elabo- 
rating on the mechanisms of geometrical frustration and 
their effects on layer configurations. In addition, we ex- 
plicitly consider the effect of the energetic cost of bending 
the layers to lie on the surface. This extrinsic bending 
effect leads to a purely geometrical ordering field for the 
layers. Though much effort has recently been directed to- 
ward this problem, particularly theoretical (2ll. [251I26L [27j 
and numerical [28l. l29l. [sol [sil . [3^ . this extrinsic curva- 
ture effect has received relatively little attention. In fact, 
bending energies can play an important role in ordering 
stripe patterns over longer length scales than would oth- 
erwise be possible [2l[ . This observation suggests the use 
of surface curvature as a motif to control self-assembly of 
block copolymers [s^. Consider, for example, the prob- 
lem of repulsive semi-flexible polymers on a cylinder. In 
the fixed density ensemble, the polymers will attempt 
to adopt an equal spacing. As shown in Fig. 1, there 
are many such structures ranging from "pin stripes" to 
"rugby stripes" , as well as the continuous class of "barber 
pole" textures. What do these have in common? They 
all have equally spaced polymers and, from the point of 
view of intrinsic geometry, they are all straight, that is, 
they are all geodesies. This can be seen, for instance, by 
cutting open the cylinder and laying it flat - all of these 
stripe textures map into straight lines on the flattened 
cylinder. However, these polymer-like lines are three- 
dimensional and their energy is a function of their em- 
bedded conformation in K^. The "pin stripe" texture in 
Fig. [T] is the only one for which the polymers have van- 
ishing three-dimensional curvature and will thus be the 
ground state. 

To this energetic accounting, we add a compression 
energy which favors equal-spacing between the lines and 
an intrinsic curvature energy within the local tangent 
plane. We write the total free energy schematically as: 

F = [compression energy] 

+ [intrinsic bending of lines] 

+ [extrinsic bending of lines] (1) 

In Sec. IIIBI we will develop explicit, geometric expres- 
sions for these three contributions. However, to frame 
our analysis, we pause to discuss the subtle frustration 
which arises from these competing terms. The compres- 
sion term measures the deviation of the layers from being 
equally spaced. On a curved substrate it is still possible 
to deflne "equal-spacing" : at each point on a given line 
we move along the geodesic tangent to the line's nor- 
mal for a fixed interlayer distance (see, for instance, the 
right side of Fig. [3]) . This construction generates equally- 
spaced layers. However, in contrast to a cylinder, on a 



curved substrate this is not straightforward. To see this, 
we consider lines on cones; the cone is flat everywhere 
but for its apex where it has a concentration of Gaus- 
sian curvature. In Fig. [5] we show three different smectic 
complexions on a cone. Each is generated by flrst trac- 
ing a smectic pattern on a flat sheet. Removing a wedge 
from the sheet allows us to construct a cone - a surface 
with vanishing Gaussian curvature everywhere except at 
the point. Geodesies on the cone become straight lines 
on the original surface and vice versa. On the left we 
see equally-spaced lines (solid) and their normal flows 
(dashed). When put onto a cone, the geodesies remain 
equally spaced. Since the intrinsic bending vanishes, by 
definition, for geodesies, this pattern has no contribution 
from the first two terms almost everywhere. However, 
due to the Gaussian curvature conentrated at the tip of 
the cone, the global effect of the Gauss-Bonnet theorem 
(discussed below) is that a sharp, intrinsic bend must be 
introduced into the layers along a seam of the cone to ac- 
count for the angle deficit. Moreover, though these lines 
are flat from the point of view of the surface, they are 
bent in and thus contribute to the extrinsic curvature 
energy as well. Since the directions along which the ex- 
trinsic curvature vanishes depends on how the surface is 
bent into three-dimensions, this effect leads to an extrin- 
sic field which can align the layers (see below). It is in- 
structive to consider the other two configurations in Fig. 
[2] as well. In the center figure a set of concentric circles 
is wrapped onto a cone. From the construction we know 
that the resulting lines are equally-spaced and have no 
sharp bends. The circles have both intrinsic and extrinsic 
curvature; we can see this by noting that the normal to 
each circle points towards its center, in the same plane 
as the circle. However, the normal to the cone's surface 
does not lie in that plane and so the circle's curvature 
has components both along the surface normal and per- 
pendicular to it. The former measures the extrinsic or 
normal curvature, the latter the intrinsic or geodesic cur- 
vature. Finally, on the right, we show a set of lines which 
all converge to a point. They are not equally-spaced on 
either the plane or the cone; by construction they are 
geodesies and for this particular embedding of the cone, 
they have vanishing extrinsic curvature. The simple cone 
geometry highlights the geometric frustration embodied 
in simultaneously minimizing the three competing terms 
in the free energy. 

In section II we introduce the geometric tools neces- 
sary to characterize stripes and develop the free energy 
that controls smectic order on curved surfaces. In sec- 
tion III we present a number of examples that illustrate 
the effects of both intrinsic and extrinsic geometry on 
striped pattern formation. Because of the many different 
limits that might be considered in this frustrated system, 
we focus primarily on zero-strain complexions where the 
lines are equally spaced (compression energy dominates) . 
This limit is closely related to the geometric optics limit 
of light propagation, as we will exploit. As we shall see, 
even with this constraint the competition between the 
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FIG. 1: From left to right: "pin stripe," "rugby stripe," and "barber pole" textures on a cylinder. These can be constructed 
from straight lines on a plane by identifying two sides of a rectangle (red). Though the layers are always geodesies, the layers 
of the left cylinder ("pinstripe") are also straight. When K = 0, the layers preferentially lie along the principle direction of 
zero curvature. 





FIG. 2: From left to right: "equally-spaced geodesies with sharp bend," "equally-spaced lines with no bend singularities," and 
"geodesies with vanishing extrinsic bending" textures on a cone. These are all constructed by taking a pattern on a flat sheet, 
as could have been done for the cylindrical texture of Fig. [1] A wedge angle is removed as shown. On the left cone the layers 
have equal spacing everywhere but for the line of dislocations along the seam and the layers have some normal curvature. In 
the center, the layers are equally spaced but there is both geodesic and normal bending. The right-hand cone has vanishing 
bending energy of all kinds, but the layer spacing is incorrect almost everywhere. The dashed lines on the cone to the left 
indicate the layer normals. 
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remaining two terms is quite subtle. In section IV we 
discuss local mechanisms (related to extrinsic curvature 
energy) which lead to long range, two-dimensional or- 
der on curved surfaces, bypassing the standard Coleman- 
Mermin- Wagner conclusions about thermal fluctuations 
in flat space [s^ . Iss! ] . Finally, we summarize our results 
and discuss open questions. 



II. LAYERED STRUCTURES ON A CURVED 
SUBSTRATE 

A. Geometrical background 

Given the focus in general relativity and string the- 
ory on intrinsic geometric concepts such as conformal 
invariance, diffeomorphism invariance, and modular in- 
variance, it is worth remembering that fascinating ma- 
terials exist in three dimensions where not only intrinsic 
but also extrinsic geometry plays a key role. With extrin- 
sic curvature in mind, recall that an embedded surface 
X(mi,U2) S R'^ has an induced metric gij (flrst funda- 
mental form), a unit surface normal N, and a curvature 
tensor (second fundamental form): 



9ij 
N 



l^iX X d2X\ 



(2) 



where the last equality follows from differentiating N • 
diX = 0. Note that our sign convention for reflects 
the standard sign for the curvature of a three-dimensional 
space curve, i.e. a circle will have a positive curvature 
when its normal is chosen to point inward. Here and 
throughout we refer the reader to [s^ for more technical 
details. 

The two principle curvatures ki and K2 are the eigen- 
values of g~^L, where the indices associated with g^j 
and Lij are suppressed to simplify the notation. The 
unit eigenvectors associated with the matrix L are the 
principal directions ei and 62. On a general sur- 
face the curvatures and directions change from point to 
point. With stripes, polymer strands, or smectic lay- 
ers in mind, consider a curve embedded in the surface, 
R(s) — X ((Ti(s), (12(3)), with s measuring the arclength. 
The curve has unit tangent T(s) — dH/ds with derivative 
dT/ds — K(s)Ncurvo(s) - an equation which defines the 
unit normal to the curve, Ncurvc(s), and the curvature, 
k(s). Because N^urvc and the surface normal N defined 
above are not necessarily at right angles, we write 



ds 



= KNcurvc = K«N + Kg (N X T) , 



(3) 



SO that the change in T is decomposed into a vector 
along the surface normal, N, and a vector in the surface, 
N X T. The coefficients k„ and Kg are the normal and 
geodesic curvature, respectively, and obey k'^ = + k'^. 



A geodesic on a surface is a curve for which the change in 
the tangent vector has no components in the surface, that 
is Kg = 0. However, even a great circle on a sphere curves 
in three-dimensional space. At any point on a curved 
surface, we can choose coordinates so that = diX are 
orthonormal so that gij 



— Sij and 



g ^ L = L ^ KiBieJ + ^26262 , 



(4) 



where ^ indicates the transpose [3(|. It follows that the 
normal curvature can be further decomposed as 



dT _ ^ dN 

ds ds 
Ki cos^ P + K2 sin^ /3 



-T • (T • V) N = T^LT 

(5) 



where T • ei = cos f3 gives the angle P between the curve 
and one of the principle directions. The last equality in 
the first line follows from the definition of L in ([2]) and 
from expanding T in terms of the two tangent vectors e; 

M- 

On a radially symmetric surface (such as the Gaus- 
sian bump discussed below) with height function h{r), 
X{r,0) — [r cos 6*, r sin 6*, /i(r)]. The principal directions 
are, by symmetry, along f and and along 9. It follows 
that if /? is the angle between the curve and the radial 
direction, then k„ = Kr cos^ P + Ke sin^ P where 



{drhy 
drh 



3/2 



K0 



r l + {drhy 



1/2 



(6) 



When KrK0 < (on the fianks of a Gaussian bump), there 
will always be an angle P for which k„ vanishes, and so 
the extrinsic geometry and extrinsic bending energy of 
stripes on a surface set a natural, local direction for the 
polymers. Now our discussion for smectic layers on cylin- 
ders (recall Fig. [T]) can be made precise. On a cylinder 
of radius — and Kg = 1/R so k„ = only when 

T • ee = 0, which selects the preferred minimum energy 
configuration. 

The cylindrical ground state is unambiguous because 
all of the textures in Fig. 1 are composed of lines drawn 
along geodesies. Note that this extrinsic effect is sensitive 
to an intrinsic quantity, the angle p. On a smooth bump, 
with a nonzero, spatially varying Gaussian curvature, the 
problem is complicated by the intrinsic geometry of the 
curves defined by the layers. Now we must also consider 
Kg, the geodesic curvature. As we will see in the follow- 
ing, the surface geometry may not only favor line con- 
figurations for which k„ 7^ but, more importantly, the 
surface geometry can prevent the stripes from achieving 
both equal spacing and vanishing Kg. 
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B. The Smectic Energy on Curved Substrates 

Lamellae, both on surfaces and in the plane, are con- 
veniently represented by the level sets of a function 
<I>(mi,M2) = *&(u) = an, where a is the equilibrium layer 
spacing and n is an integer labeling the layer. The phase 
field <& is related to the average mass density p(u) of the 
layers by p(u) = + pi cos(27r$/a). The gradient 9^$ 
contains information both about the layer spacing, re- 
lated to its magnitude, and the direction normal to the 
layers. We will use V to denote the covariant deriva- 
tive on the surface, defined as Vi<i> — 9^$ for a scalar 
function <I>, and V^u^ = diV^ + F^^.w'^ for (contravari- 
ant) vectors [s^ . When unambiguous, we will also use 
• to indicate contraction: thus, v • V$ = v^di^. Note 
that V • V$ gives the change in $(?/) when one moves in 
an arbitrary direction v on the surface. Upon specializ- 
ing to the case where v = V$, the resulting change is 
jV'l'p — V$ • V$ — g^^di^dj^. Hence, equal layer sep- 
aration implies the condition |V$p = 1. The strain, e, 
measures the deviation from equal spacing. While there 
are many possible forms for e, depending on microscopic 
details, we only require that e vanish when |V<i>p = 1 
and, in the small deformation limit where the Eulerian 
displacement is m = z — $, that e ^ dzU where z is the 
local direction of the layer normal. A suitable form is 



1 



|V$|' 



(7) 



The energy cost of small deformations from equal spac- 
ing will be an expansion in powers of e, with the lowest 
order term being e^. If the layers are equally spaced, 
then 1 = |V$p and, upon differentiation, we find that 
= (V$ • V) V$ = (n • V) n where = g'^ dj<l>/\V^\ 
is normal (within the tangent plane) to the curve used 
to define a layer. The condition (n • V) n = is pre- 
cisely the geodesic equation for the layer normal. We see 
then, independent of the form of the strain, when e = 
(thus, minimizing the strain energy) the layers are spaced 
evenly along geodesic curves. If we consider only zero 
strain solutions, this condition reduces the problem of 
finding smectic textures to a first- order evolution equa- 
tion. Thus, the normals lie on "straight lines" on the 
surface {i.e. geodesies); layer curvature leads to their 
convergence into singularities. 

Indeed, if ^ measures the distance between two 
geodesies along a third, perpendicular geodesic, the best 
local approximation to parallel trajectories leads to (36| 



ds2 



-7^ = 



(8) 



where K is the Gaussian curvature. This geodesic devi- 
ation equation naturally reminds us of geometric optics 
[sj , in which the layer normals act as rays and the Gaus- 
sian curvature, K{ui, U2) as a variable index of refraction. 
In regions of negative curvature the normals diverge as 
through a diverging lens; in regions of positive curvature 



they converge. Converging patterns of rays form caustics; 
in the language of ordered lamella, caustics are curvature 
singularities in the lines. On a curved surface, a set of 
geodesies normals initially perpendicular to some layer 
may cross a finite distance away, leading to a cusp-like 
boundary where the geodesies finally converge. This sys- 
tem thus provides a low-dimensional analog of gravita- 
tional lensing ^38|, |39i] . Cusps are regions of high bending 
energy and thus, in the true ground state, they will be 
smoothed out at the cost of compression energy. 

We are now in a position to formulate the complete 
free energy: 



F 



B 
"2 



dA 



(9) 



where B is the bulk modulus and the couplings Xg and 
An are length scales, typically on the order of a column 
diameter, which measure the relative importance of bend- 
ing and compression energies. In principle Ag ^ A„: the 
A„ — > limit is the "generally covariant" limit where 
extrinsic effects are irrelevant. On the other hand, if 
we were trying to decorate the surface with semi-flexible 
polymers with persistence length Lp, we would have 



A^ 



= ksTLph/ {Ba ), where h is the thickness 



of the polymer layer normal to the surface and a is the 
average spacing between the polymers [40^. In general, 
nonzero Xg and A„ will lead to corrections to the simple 
picture of caustic singularities sketched above. 

To better understand Eq. ([9|), consider the three dif- 
ferent decorations of the cone shown in Fig. 2, assuming 
finite cones of radius R. On the left, while there is some 
extrinsic layer bending leading to an energy on the or- 
der of BX^ Ini?, the dominant energetic cost arises from 
the bend wall on the seam: there the equal-spacing con- 
dition will breakdown in order to smooth out the sharp 
kink. Equivalently, we can view this as a row of disloca- 
tions as arises in low-angle grain boundaries with elastic 
core energy per unit length scaling as Ba^. In either 
case, we see that the energy of the grain seam scales as 
Fieft ^ BaR. Upon focussing on the central cone, we see 
that the smectic complexion has energetic contributions 
from both bending moduli and Fcenter -S(Ag-l-A^) Ini?. 
Consider finally the cone on the right; the texture there 
has no bending energy at all, but a very large has com- 
pression energy, Fj-ight ^ BR^. From this rough analysis, 
we would expect that away from regions of large Gaus- 
sian curvature (in this case the conical singularity), the 
"bulls-eye" pattern of the center cone of Fig. 2 will be 
the dominant line texture. However, depending on the 
various elastic constants this can break down for finite 
regions around the the peak. In the following we will 
consider these effects, in particular in Section IV. 

For further insight into Eq. it is instructive to con- 
sider this free energy in the context of weakly deformed 
two-dimensional smectics coating a cylindrical surface 
with a flat metric and no singularitiy. We introduce the 
usual one-dimensional displacement fleld w(ui,U2) and 
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consider the simple case of a cylinder discussed in the 
introduction. In the coordinate system of Fig. [31 with 
(ui, U2) = (x, z) = r, we take as our level set function 



<i>(x, z) = z • r — w{r) 



(10) 



Consider the limit of gentle, slowly-varying undulations 
superimposed on the texture shown at the bottom left of 
Fig. [21 For a cylinder of radius R, it is straightforward 
to show that the free energy then becomes 

F=jJdA [{d^wf + Xl {dlwf + [Xn/Rf sin^ {d^ 

(11) 

where dxw{T) « /3 is the local tilt angle that the stripes 
make with respect to the x-axis. The last term in Eq. 
(fTTjl follows from Eq. ([5]), where we choose axes of prin- 
ciple curvature such that ki = and K2 = In addi- 
tion to the usual terms describing a rotationally-invariant 
2d smectic 41|, there is now a weak ordering field pro- 
portional to sin* j3. The field is weak in two senses - its 
strength is proportional to \/R^ and hence vanishes in 
the limit of a cylinder of infinite radius. Furthermore, 
it goes like j3'^ for small tilt angle deviations from the 
preferred direction. Nevertheless, this field is enough to 
break the symmetry and give rise to a preferred direc- 
tion for the lines. As discussed above, this field arises 
from the extrinsic curvature tensor and will vary spa- 
tially in both magnitude and direction on more general 
curved surfaces. In the remainder of this paper, we first 
neglect this weak field (for a Gaussian bump of height /lo 
and size i?o, h/R^ replaces 1/i?) and focus on satisfying 
the constraint of equal layer spacing imposed by the first 
term of Eq. with a simplifying boundary condition 
at infinity. We then discuss effects due to this extrinsic 
ordering term with free boundary conditions at infinity. 

An anisotropic membrane would have two bending 
moduli, one for bends along a "hard" direction, the other 
for bends along the "soft" direction. If the anisotropy had 
a vector character, we could also include a cross term, al- 
lowed by symmetry. Defining t as the unit vector lying 
along the layers, an alternate form for the normal curva- 
ture is K„ = t^PLij, where Lij is the surface curvature 
tensor given in equation ([2]). The other two allowed terms 
are n^n^ Lij and n^PLij, respectively. The magnitude 
of the moduli for these additional terms depend on the 
molecular details of the smectic layers. These elastic con- 
stants are different, in principle, from the intrinsic bend- 
ing modulus BAg in For concreteness, we will focus 
our discussion on a monolayer of block copolymer cylin- 
ders lying on the surface. The columnar phase of neat 
{i.e. monodisperse with no solvent) block copolymers, 
though bearing resemblance to the columnar hexagonal 
phase of liquid crystalline polymers, is actually an in- 
compressible three-dimensional elastic medium. Strong- 
segregation calculations (4§| suggest that the columns can 
be viewed as semi-flexible rods. In that case "bending 
along the columns" is the "hard" direction, being more 
energetically costly than "bending perpendicular to the 



columns" for a few-layer coating of a curved substrate. 
Roughly speaking, the columnar phase is similar to a 
corrugated sheet with the columns corresponding to the 
corrugations. However, as the diblock film grows, bend- 
ing along the columns leads to deformations that are in- 
dependent of the thickness since lamella can slide past 
each other with no cost. On the other hand, bending 
perpendicular to the columns requires large amounts of 
crystalline strain. For a thin film composed of only a few 
layers, the bending energy we consider here should dom- 
inate. However, how the introduction of low-angle grain 
boundaries and surface energies affects these calculations 
is an open question ^42.]. To keep our analysis from be- 
coming highly ramified, the only non-vanishing extrinsic 
elastic modulus will be bending along the column tan- 
gents. 



C. Mechanisms of geometric frustration 

With the spate of recent work on both crystalline order 
and nematic order on curved substrates, it is valuable to 
compare and contrast these systems with the smectic - 
a phase with one-dimensional crystalline order that lives 
between the crystal and the nematic. We shall see that 
the smectic presents issues all its own and affords a fresh 
arena for the interplay between geometry and soft mate- 
rials. 

A nematic liquid crystal on a surface is described by a 
unit vector n which lives in the tangent plane of the sur- 
face. The standard three-dimensional Frank free energy 

m 

F[n] = IJd^"^ {^1 • + [n • (V X n)]^ 

+i^3[(n- V)n]'} (12) 

is modified in two ways: First, n • (V x n) = when n 
lies on a smooth, two-dimensional surface and depends 
only on the two surface coordinates u = {^1,1*2}. The 
two-dimensional nematic free energy then reads 



Fnem — 



j dA{V -nf + ^ J dA [(n • V) n]' 



(13) 

where dA — d^u^ and g is the determinant of the metric 
tensor g^j. The first term in Eq. (|13p penalizes director 
splay, whereas the second term penalizes director bend. 
When Ki = = Kp the free energy in Eq. is 
isotropic and can be cast in a form that naturally lends 
geometric insight. 

Consider a local angle field /3(u), corresponding to the 
angle between n(u) and an arbitrary, orthonormal local 
reference frame whose basis vectors we label ei(u) with 
i — 1,2. In the one Frank constant approximation, the 
free energy in Eq. (|13p can be recast in the form 



dAg'^{d,e-n,){dj(3-n,) , (m) 
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where f2i(u) is a connection that plays the role of the 
ChristofFel symbols by compensating for the rotation of 
the 2D basis vectors ei(u) in direction i and W x Q = K 
[4^ . Since the curl of f2,i(u) is equal to the Gaussian 
curvature K{u), the nematic energy cannot vanish on 
a surface with non-zero Gaussian curvature. Note that 
rii(u) is, in general, a non-conservative field, so we cannot 
minimize (jl4p by setting di(3 equal to f2i everywhere on 
the surface. This property is a manifestation of a more 
general mechanism, commonly referred to as geometric 
frustration, to indicate situations where the molecular 
arrangement favored by local interactions cannot be ex- 
tended globally. 

As the Gaussian curvature of the substrate increases, 
defects are generated in the ground state to lower the 
energy cost of geometric frustration. Their energetics 
is analogous to Coulomb particles interacting with a 
smeared out electrostatic charge given by the Gaussian 
curvature. This nontrivial result can be rationalized by 
examining the free energy in Eq. ()14|1 and noticing that 
the connection r2i(u) and the Gaussian curvature K(u) 
are analogous (in two dimensional electromagnetism) to 
a frozen vector potential and the magnetic field respec- 
tively. The topological defects, i.e. disclinations, behave 
as monopoles in the dynamical field (3{u) whose inter- 
action with the geometry of the surface is mediated by 
the geometric gauge field rii(u). An additional couphng 
between defects and the metric arises from the metric 
factors that appear, for example, in the surface element 
dS = (Puy/g independently of the connection [45]. 

The physics of geometric frustration is at work also in 
the more complicated setting of curved space crystallog- 
raphy 24i, i46,]- Now the orientational order of the 
nematic is supplemented by translational degrees of free- 
dom. For gently deformed surfaces, the crystalline energy 
can be expressed in terms of the Lame coefficients /x and 
A m 




where x = {x, y} are cartesian coordinates in the plane 
and Uij{x) = ^ [diUj{x) -t- djUi{x) -I- Aij{x)\ is the strain 
tensor. Compared to its flat space counterpart, the strain 
tensor has an additional term Aij {x) = dih{x)djh(x) that 
couples the gradient of the displacement field Ui{x) to 
the gradient of the surface height function h{x). The 
field Aij (x) is a tensor version of the connection fli in- 
troduced above to describe orientational order on curved 
surfaces. Indeed the curl of the tensor field Aij{x) is 
equal to the Gaussian curvature of the surface K{x) — 
—£ii£jkdidkdih(x)djh{x) where e^j is the antisymmetric 
unit tensor {cxy = —Syx = 1) [is, 49]. By the same rea- 
soning as before, the integrand of Eq. p5|) and hence 
the ground state energy cannot be made to vanish. This 
is the mathematical mechanism by which geometric frus- 
tration enters the physics of curved crystals. It can be 
grasped more intuitively by recalling that bending a plate 



/ 




FIG. 3: Left: Turning angles A9i defined on a triangle. 
Right: A square with analogous turning angles formed from 
two uniformly-spaced layers (solid) and two normals (dashed). 
The normals are geodesies so Kg = along them. 

into a surface of non vanishing Gaussian curvature nec- 
essarily causes it to stretch p7| . 

Smectic liquid crystals, on the other hand, can main- 
tain uniform layer spacing (and hence achieve zero strain) 
even in the presence of Gaussian curvature. The Gaus- 
sian curvature, nevertheless, couples to the curvature of 
the layers. To see how this coupling appears, we start 
by noticing that the bend coupling proportional to K3 in 
equation (fT3| resembles the geodesic equation for curves 
tangent to the director. We can establish an intuition for 
smectic patterns, therefore, by studying ncmatics with 
very large if 3, and identifying the nematic order parame- 
ter as the layer normal. Nematics on spheres are required 
to have a net +2 topological charge, which tend to break 
up into four -1-1/2 disclinations. When K3 = Ki, the 
disclinations lie on the corners of a tetrahedron [H^, |5l[ ; 
when K3 is large, however, the -1-1/2 disclinations lie on 
an equator [s^] and the local texture is a lines of longi- 
tude structure. At the other extreme, when Ki is large 
the director takes on a lines of latitude texture. 

Since there is no difficulty to finding geodesies in any 
direction at any point on a curved surface, there is no 
local obstruction to constructing equally-spaced layers 
on any surface. However, doing so while simultaneously 
finding layers which are also geodesies is impossible. This 
can be seen by constructing a rectangle with two oppo- 
site sides given by adjacent layers and the remaining two 
sides given by geodesies normal to both layers (see Fig. 
[3]) . Upon applying the Gauss-Bonnet theorem [36] to this 
contour, we obtain 

J dAK = 2tt -Y^AOi- <j) dsKg, (16) 

where A9i = tt/2 is the turning angle at each corner of 
the rectangle, we find 

J dAK = J ds Kg ~ J^ds Kg. (17) 

The integral on the left-hand side is over the area of the 
square while the integrals on the right-hand side are over 
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the two adjacent layers. Since / dAK 7^ in general, it 
is not possible for = on both layers. 

Despite our local ability to set e = in equations ([7|) 
and there may be global obstructions to finding low 
energy configurations with vanishing compression strain. 
To see this, we first recast ([T6|) locally. Upon using the lo- 
cal coordinates {x,y) and introducing the surface metric 
gij, Eq. pT|) becomes 



dxdyy/gK = 



dn 



dn I ds Kn 



(18) 



where d/dn is the derivative along the normal direction 
and we have used the fundamental theorem of calculus. 
After rewriting the integral on the right in the local co- 
ordinates, we come to the local relation 



(19) 



Recall that the geodesic curvature is the fractional rate of 
change of the length of an arc as it is moved perpendicular 
to itself and that positive curvature implies that normal 
evolution shrinks the curve [s^. It follows that Kg = 
—g^^^'^dny/g where we define the normal direction to be 
along the curve's (inward) pointing normal. 

Using Eq. we find ^^K = {dny/g)Kg + y^dnHg. 

Dividing through by y^, we find 



dnKn 



K. 



(20) 



This evolution equation encapsulates the geometric frus- 
tration implicit in the Gauss-Bonnet theorem. Indeed, 
following the same reasoning, Eq. ([8]) is also a local form 
of Gauss-Bonnet; we may take an arc of length ^ connect- 
ing the two diverging (converging) geodesies at s. Then 
the curvature of those arcs is Kg = —^~^{d^/ds), and it 
follows that —K^ = d^S^/ds^ — —dKg/ds£, — Kg{d£,/ds), 
which is identical to ((20|l . 



III. SMECTIC SCATTERING FROM CURVED 
SURFACES 

In this section, we review and generalize our prior re- 
sults on smectic textures on a simple Gaussian bump [2l| 
to more complex geometries. 



A. The Gaussian Bump 

Consider smectic order on a curved substrate described 
in the Monge representation by a height function h{x,y); 
the Gaussian curvature is then given by ^] 



K 



dlhdlh - {d^dyhf 



1 + [d^hf + {dyhy 



(21) 



If the whole surface can be described by h{x,y) then 
it is topologically equivalent to the plane and, by the 



Gauss-Bonnet theorem, / dAK = 0. We start with the 
Gaussian bump used to study nematic order in Ref. [5l| , 
h{r) — hoe^^ ^'-^^ \ with = x^ + y"^. K is positive 
near r = and negative for large r. At an intermediate 
radius r — R, K vanishes. A particularly simple smectic 
configuration results from choosing radial geodesies as 
the layer normals, in which case we would find uniformly- 
spaced, azimuthal layers, a disclination at the top of the 
bump, and a power law decay of Kg away from the center. 
Is it possible, however, to generate a configuration free 
of topological defects on a bump? 

Suppose we start with straight layers at a; = —00, so 
that the layer function ^{x, y) defined in Sec. IIIBI obevs 
^ X as X —^ —00. This boundary condition would 
describe an experimental setup where the layers of a di- 
block columnar phase grow along a temperature gradient 
parallel to x. We expect layer by layer growth nucleated 
from a boundary at large negative x. The geometry of 
the substrate leads to the formation of singularities, but 
rather different than the isolated disclination with con- 
centric circular layers centered on the bump discussed 
above (see Fig. |4]). We call these singularities caustics 
because, just as in geometrical optics, these are places 
where many initially parallel light rays or, in our case, 
geodesies converge. At these locations the value of $ is 
well-defined but V$ is discontinuous. Though it might 
be tempting to call these cusps "defects", they are not. 
Recall that near a dislocation the smectic order param- 
eter vanishes and the phase field $ takes on all values 
around the defect. Similarly, in the vicinity of a discli- 
nation, the nematic order vanishes and the layer normal 
takes on all directions at the defect. In contrast, the 
cusp singularities in Fig. |4]have definite values of $ and 
V$ does not wind through all possible directions around 
these singularities. Thus, although there is a disconti- 
nuity in V<I>, it is not of the same nature as that of a 
disclination. In the analogy to optics, we would say that 
a dislocation is a place of vanishing amplitude, while a 
caustic is a location of very high, if not infinite, ampli- 
tude [1^. Figure m shows a birds-eye view of a Gaussian 
bump coated with stripes, with the geodesies as dashed 
lines and the layers themselves as solid lines. The red 
circle indicates the locus of points for which K = and 
we have shaded the regions according to the magnitude 
of Kg . We could also have predicted these cusps by deter- 
mining where the curvature of the smectic layers diverges 
along the geodesies defined by the layer normal. To do 
this, we integrate (PD|) along the normal geodesic pass- 
ing over the top of the bump. We find that Kg diverges 
a finite distance past the center of the bump, shown in 
Fig. m this divergence indicates the onset of the infinite 
curvature cusps in the layer lines. 

From a more global perspective, the cusp angle is a 
measure of the integrated Gaussian curvature. To see 
this, we again use the Gauss-Bonnet theorem to describe 
evolution of the cusp angle as a function of distance from 
the center of the bump. The key is to define a geodesic 
triangle with one edge along the x-axis, another edge 



FIG. 4: Layers (solid lines), normals (dashed lines) and 
geodesic curvature (color) for a bump with aspect ratio 
ho/R = 3 projected onto the xy— plane. The circle of zero 
Gaussian curvature is depicted in red. A scale bar of length 
R has been provided. With the constraint of equally-spaced 
layers along the normals, the curvature induces the formation 
of a grain boundary (heavy dashed line) that extends infinitely 
far to the right of the bump center. The apparent unequal 
layer spacings in the figure are an artifact of the projection. 
The inset shows the bump, in perspective, overlaid with a 
square grid. Note that the grid lines are not equally-spaced 
in this case. 



along the layer parallel to the y-axis at x — s- —oo and 
a third edge along a normal curve, which is a geodesic 
by construction. The resulting geodesic triangle has two 
exterior (or interior) tt/2 angles. Denote the remaining 
internal angle by a. The corresponding external angle, 
IT— a, is constrained by the Gauss-Bonnet theorem, which 
leads to 



dAK = 2tt- 



2 + 2+(--") 



(22) 



where the integral is taken over the inside of the trian- 
gle. Notice that, as the cusp location x becomes large 
and positive (see Fig. [4]), a{x) since the triangle 
incorporates an increasing amount of area in the entire 
upper half-plane. Because the integrated Gaussian cur- 
vature is zero in the upper half-plane, the angle a nec- 
essarily decreases with increasing x. Though the grain 
boundary persists infinitely far to the right of the center 
of the bump, the cusp angle asymptotically vanishes. In 
the case of a Gaussian bump, the angle will fall off as 
Q-(r/R) ^ though the details will vary for other surfaces. 
Similar to geometric optics, the formation of these cusps 
should not depend sensitively on the exact geometry of 
the substrate; rather it is a function of the topology as 
characterized by the intrinsic curvature. 



B. Smectic optometry 

We have considered other substrates, all topologically 
equivalent to the plane so that again, the integrated 



FIG. 5: A converging lens: Layers (solid lines, only 
upper half-plane shown) and normals (dashed lines, only 
lower half-plane shown) for a bump with h = ho[{x^ + 
y^)/i?]exp[-(r/i?)V2) {ho/R = 1). A scale bar of length 
R has been provided. All quantities exhibit mirror symmetry 
about the horizontal midline. Boundaries between regions of 
Gaussian curvature with different signs are delineated in red. 
Grain boundaries, the analogues of caustics in geometrical op- 
tics, are shown as bold dashed lines. Note the focusing of the 
normal lines onto the grain boundary in the lower half-plane. 
All lengths are measured in units of the bump width R. The 
surface is shown in perspective in the inset. 



Gaussian curvature vanishes, / K dA — 0. In Fig. [5] 
we depict both the geodesies and layers generated from 
the same boundary condition as a; ^ — oo that we consid- 
ered in the last section; now, however, we have chosen a 
more complex axisymmetric bump with multiple regions 
of positive and negative K. Note that geodesies that do 
not go through the central region of the bump, outlined 
by the outermost red circle, will generate a set of cusps 
at large x along the midline similar to those in Fig. H) 
However, we now find that there are additional caustic 
lines that form as a result of the focussing from the in- 
verted dimple of the bump. These "fold" caustic lines, 
in this case, converge toward the x-axis and eventually 
end there. To see this, we adapt the discussion in the 
last section: geodesies that avoid the central region en- 
tirely lead to larger-angle cusps on the x-axis than those 
that go through the same region. At large x, the pattern 
must be that of the simple bump in Fig. |4]and thus there 
will be no auxiliary caustics at large x and y. It follows 
that any extra caustics must converge to the x-axis. We 
can see this in another way by considering the focussing 
equation ([8]). Pairs of geodesic normal lines that remain 
in the outer region always diverge, because K < 15^. 
However, a pair of geodesies, one of which remains in the 
outer region and which enters the inner region will still 
diverge less slowly and, indeed, pairs which both travel 
through the annulus with K > will converge. Hemmed 
in by the pairs that are always diverging, we see that 
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FIG. 6; Layers (solid lines, only upper half-plane shown) 
and normals (dashed lines, only lower half-plane shown) for a 
saddle-bump with h{x, y) = 0.6[(a;^ - y'^)/R] exp[-(r/_R)^/2]. 
A scale bar of length R has been provided. Boundaries be- 
tween regions of Gaussian curvature with different signs are 
delineated in red. Grain boundaries are shown as heavy 
dashed lines. The surface is shown in perspective in the inset. 



these additional cusps must form and that the "outer" 
geodesies overtake the "inner" geodesies, leading to two 
extra cusp line grain boundaries, one of which is shown 
in the lower half of Fig. 

Just as microscopists design and choose multiple lenses 
to manipulate an image, one might hope to engineer 
a surface whose geometry leads specific smectic tex- 
tures. To further our intuition, we have studied the non- 
axisymmetric saddle-bump shown in Fig. [SI with the 
polynomial — multiplying a Gaussian envelope. As 
shown in the inset, this surface has a saddle in the cen- 
ter, surrounded by four lobes of positive curvature. Upon 
choosing the same boundary conditions as x —^ —oo as 
before, we find a line of cusps along the x-axis and two ex- 
tra off-axis caustic lines that disappear as they approach 
y = 0. In order to find an example for which the auxiliary 
lines diverge (a diverging instead of converging lens), we 
consider, for instance, the same saddle-bump structure, 
but rotated 45° so that the geodesies that start near the 
X-axis minimize their transit across regions of positive 
curvature. These geodesies are always diverging, picking 
up only a small amount of ii' > as they exit the central 
saddle region. As a result, there are no cusps on the x- 
axis - the layers are perfectly smooth since the geodesies 
associated with the layer normals are not focused to this 
line. As shown in Fig. [71 one can see that the cusps in the 
layers form symmetrically, off-axis and that, for the bump 
shown, they bend away from the line y — 0; it would be 
impossible for them to converge since the layers are reg- 
ular on the X-axis. It seems likely that, by adjusting 
the saddle-bump parameters, the lines of cusps could be 
made parallel to x. These simple, high-symmetry cases 




FIG. 7: A diverging lens: Layers (solid lines, only up- 
per half-plane shown) and normals (dashed lines, only 
lower half-plane shown) for a bump with h{x, y) = 
0.55{xy / R) exp[— (r/i?)^/2]. A scale bar of length R has been 
provided. Boundaries between regions of Gaussian curvature 
with different signs are delineated in red. Grain boundaries 
are shown as heavy dashed lines. The surface itself is shown 
in the inset. 

provide intuition for how layer focussing could be used 
to construct desired patterns far from the bump. It is 
amusing to consider other boundary conditions (for in- 
stance, a circular or parabolic arc at infinity incident on a 
bump with little or no symmetry. Certainly the "optical" 
elements that we have studied here could, themselves be 
used to control the layers and geodesic normals if they are 
put far enough apart; the bump in Fig. [Sjacts as a con- 
verging lens, while the bump in Fig. [Tjacts as a diverging 
lens. The effect of such a compound "lens" can be seen 
in Fig. [H The converging bump on the left leads to the 
formation of an x-axis, downstream grain boundary, as 
it does in Fig. [51 this particular the grain boundary ends 
due to the diverging effect in the negatively-curved center 
region of the bump from Fig. [71 but two new diverging 
grain boundaries appear downstream of this element. We 
look forward to a future where "smectic optometrists" 
manipulate the layered order with Gaussian curvature to 
make novel devices and materials. 



IV. A LOCAL ORDERING FIELD FROM THE 
EXTRINSIC GEOMETRY 

A. Ground States: Involutes and Evolutes 

In currently available experimental data [l^, HH , the 
stripes form all at once across the surface as the block 
copolymers cure. Similarly, if a nematic on a sub- 
strate were cooled into the smectic phase, the forma- 
tion of striped order might not proceed from left to right 
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FIG. 8: Smectic optometry: a compound "lens" built from the bumps described in Fig. [5] and [7] Layers are again denoted 
by solid lines and layer normals by (light) dashed lines. A scale bar of length R has been provided. The grain boundary that 
appears due to the convergence of the left-hand bump ends at the diverging bump on the right. A second pair of caustics/grain 
boundaries appear to the right of this bump. Only one is shown, as the heavy dashed line. 



as imagined above unless a temperature gradient or a 
strongly nucleating boundary condition at x — > — oo were 
imposed. Moreover, in the prior subsections, we have 
only considered the effect of the intrinsic geometry, em- 
bodied in the Gaussian curvature, on the formation and 
orientation of the layered structure. This approximation 
amounts to taking the limit A„ and Ag to zero in Eq. 
We know from our discussion in the introduction that, 
when the bending couplings A„ and Ag in Eq. ^ are 
important, the layers prefer to be perpendicular to the 
lines of if = in Figs. |4][7l Can a smectic texture with 
equally-spaced layers be found so that it agrees with the 
preferred direction along lines given by the K = locus? 
In this section we will consider the additional effect of 
the normal curvature coupling A„ but set Ag = 0. The 
intrinsic curvature term differs essentially from the other 
two in ([1]); equal spacing is a local constraint and the ex- 
trinsic curvature generates a local preferred direction as 
illustrated in Eq. pT|) . However, the intrinsic or geodesic 
curvature can always be set to zero locally because there 
is a geodesic pointing in every direction at each point 
of the surface. As illustrated in Fig. [3] it is only after 
evolution normal to the stripes that the geodesic curva- 
ture builds up from zero. Thus, the geodesic curvature 
term does not set a local directional contstraint. For this 
reason we neglect Ag in the following discussion and con- 
cern ourselves with the orienting effect arising from A„. 
Presumably, local ordering fields would dominate the ki- 
netics, particularly in layer-by-layer growth. 

For simplicity we return to the simple Gaussian bump 
in Fig. m First, consider the region for which K < 0. 
Recall from the discussion around Eq. ^ that there 



always exists a critical angle f3 relative to the principal 
directions such that k„ in Eq. ^ vanishes. 

As discussed in Sees. I and II, the layers will point 
in the radial direction along the ring where A' = to 
minimize the normal bending energy. Remarkably, it is 
possible to construct an equally-spaced array of lines sat- 
isfying this boundary condition with free boundary con- 
ditions at infinity. This texture is given by the involutes 
of the ring of zero Gaussian curvature, generalized to a 
curved surface. Recall that in flat space, the involutes of 
a closed plane figure can be generated by wrapping an 
inextensible string many times around the figure, attach- 
ing a pen to the end, and then unwrapping the string 
under tension ^3d\. In optics, the evolutes are the wave 
fronts which form a caustic ~ the patterns which we have 
shown in Figs. |4][7] are the evolutes, generalized to a 
curved surface, of the straight line at x = — oo. The gen- 
eralization is straightforward: we now require that the 
string lie within the surface so that distances are mea- 
sured intrinsically. 

We can make this observation more precise by find- 
ing a family of geodesies normal to the layers, from 
which we can generate uniformly-spaced stripes. On sur- 
faces of revolution, these geodesies are characterized by 
Clairaut's theorem |36|, rs'mO — tq for a constant tq, 
where is the angle the geodesic makes with respect to 
the radial direction. The layers, which lie perpendicular 
to the normals, therefore make an angle — tt/2 — 9 with 
the radial direction so cosip — ro/r. As for the bumps 
discussed above, r is the distance from the axis of revolu- 
tion; that is, the radial projected distance. When r — tq, 
corresponding to the circle if = of the Gaussian bump. 
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FIG. 9: This texture of involutes generalized to a curved sur- 
face has both uniform layer spacing and minimizes the normal 
curvature along the circle 7^ = 0. Note that the layers bend 
tightly to be perpendicular to this circle, leading to a diverg- 
ing curvature. To the right is a top view: a projection onto 
the {x, y) plane. 



-0 = and the layers are radial; when r ^ oo, — tt/2 
and the layers are azimuthal. Thus, in the absence of 
defects, the extrinsic alignment forced at if = leads to 
a circular "lines of latitude structure" for the layers as 
r — !■ cx). In Fig. [5] we show these involutes drawn on the 
Gaussian bump of Fig. We also note that just as the 
layer curvature diverges on the caustic cusps in Figs. [4][8l 
the curvature of the layers diverges on the ring. To see 
this directly, calculate 



1 



Kg = - V • n = - —^di [^fgn'] 
v9 



(23) 



where the layer normal has components n'' — 

simp/ \Ji + (drh)^ and n*^ = cosip/r. For the layers de- 
fined above, this yields. 



(24) 



which diverges as r 



\ Note that our involute ansatz 



interpolates between the minimizer of normal curvature 
near K = and the low energy, azimuthal "bulls-eye" 
pattern of the central cone in Fig. [2] Though it might 
be natural to assume that the bulls-eye pattern would 
minimize the overall energy, it is frustrated by the ori- 
enting normal curvature energy at the ring of K = - 
experiment also shows that the layers lie perpendicular 
to this ring [ll[2l|. 

Certainly in the case of the involutes, the diverging 
curvature at if = will be softened by breaking the zero 
compressional strain condition. From dimensional anal- 
ysis one can argue that at length scales longer than A„ 
and Ag, the strain energy dominates the curvature energy 
and that the e = approximation is valid. However, at 
these cusp-like singularities, the curvature always domi- 
nates the energy no matter the scale of A„ and Ag. The 
singularities at these cusps are relaxed by the competi- 
tion between compression and curvature elasticity. The 
proper analysis of the breakdown of equal spacing re- 
quires a more complete theory which also allows the for- 
mation of edge dislocations in the smectic order [s^ . 



B. Local Order 

Equation (O, when inserted into Eq. ([9]), endows the 
layers with a preferred directionality by minimizing k^, 
and, hence, the normal curvature energy. In regions 
of positive curvature, the minimal energy direction lies 
along the principal direction with the smallest curvature. 
In regions of negative Gaussian curvature, ki/k2 < 0, so 
we can solve equation ([5]) for k„ = 0, finding a preferred 
angle f3o given by (see Eq. (O) 



tan/3o = \l -■ 

K2 



(25) 



It is instructive to again consider the case of a cylinder 
of radius R, in which ki = along the cylinder axis 
and K2 — l/R. The preferred angle is /3 = which im- 
plies that the layers all lie along the cylinder axis. Why 
this should be so is immediately apparent from Figure 
[H when the layers lie along the cylinder axis they are 
straight in three dimensions whereas when they are az- 
imuthal, each layer has curvature 1/i?. 

On our Gaussian bump, h{r) — hoe"^ and there 

is a ring of if = 0; near this ring, the surface is cylinder- 
like and we should expect a preferred radial direction for 
the layers as discussed above. Outside the ring, K < 
0, and the preferred direction follows equation in a 
surface dependent way. There is no general behavior that 
can be inferred via the expressions for the curvatures in 
Eq. ^ when if < - the precise form of h{r) dictates 
the ordering direction. Even far from the bump there is 
no universal direction; any angle is possible depending 
on lim^^cx) Kr{r)/Ke{r). 

To probe whether these low curvature directions are 
actually adopted by the layers, we compute the normal 
bending energy cost for deviation from the preferred di- 
rection for a small patch of layers from Fig. [9] in the re- 
gion where if < 0. A straightforward calculation yields. 



S£ 



(Kg - Kr) 

-BXlK 



dA sin^ (2/3o) {Sf3f 
dA {5l3f 



(26) 



when /3o > 0. We have used the condition k„ = and 
inserted the preferred angle (3q in order to simplify Eq. 
(I26p : the second equality follows from trigonometric iden- 
tities. The magnitude — if > decreases as r ^ oo but 
is order up to r/R k, 3. At small angles, however, 
sin(2/3o) is small and the ordering field is particularly 
weak. Near regions where /3o = it follows that = 
and 



5E 



BXt 



4 / dA {Spy 



(27) 



Thus there is still a weak, anharmonic ordering field even 
as the layers approach the ring of zero Gaussian curva- 
ture, just as for the cylinder described by Eq. pT|) . 
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V. CONCLUSION AND OUTLOOK 

We have outlined a theory of uniformly-spaced lines 
on a curved surface and used this to understand smectic 
liquid crystal textures on curved substrates. We started 
with the more general free energy of Eq. ([9]), but special- 
ized to the case where the compression energy dominated 
(Ag = \n = 0). We then studied the effect of the weak 
ordering field embodied in a small A„ ^ Ag. This the- 
ory is geometrical in nature, and demonstrates that sub- 
strate Gaussian curvature subtly frustrates smectic order 
by bending layers and leading to defects and grain bound- 
aries. We have neglected kinetic effects [53| focussing on 
ground state configurations. 

We now comment on additional limitations of our 
approach for real experimental systems. In particular, 
two-dimensional smectic order is destroyed by thermal 
fluctuations and the nucleation of defects The 
remaining intermediate-range order leaves a labyrinthine 
phase of stripes with anisotropic correlation lengths ^| 
and fx that can be considerably larger than a layer 
spacing, as occurs on a flat substrate. Thus, intrinsic 
curvature effects (the only effect possible in a flat 
system) can be weak, and are likely to be less important 
than extrinsic curvature in determining layer structure. 
When the effects of normal curvature are compatible 
with uniform spacing, we expect the layers to be very 
well-ordered over long length scales precisely because 
the curvature acts as an "ordering" field. This induced 
intermediate-range order on a corrugated surface does 
indeed seem evident in experiments of block copolymer 
cylinders on a bump p^ . [2l| . In particular, the layers 



will transition from being radial near the radius of 
K = Q to azimuthal far from the bump. Where bias 
introduced by normal curvature disagrees with the 
constraints required by uniform layer spacing, we expect 
there will be defects in the layers that accommodate 
the extrinsic curvature effects more closely. Therefore, a 
more quantitative understanding of the role of defects 
is necessary to build a complete picture of the ground 
state of smectic liquid crystals on curved surfaces [ssj . 
Our results can serve as a minimal template for smectic 
textures from which a more detailed theory of defects 
can emerge. 
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